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Abstract. We examine domino tilings of rectangular boards, which are in natural bijection with 
perfect matchings of grid graphs. This leads to the study of their associated perfect matching 
polytopes, and we present some of their properties, in particular, when these polytopes are Goren- 
stein. We also introduce the notion of domino stackings and present some results and several open 
questions. Our techniques use results from graph theory, polyhedral geometry, and enumerative 
combinatorics. 



1. Introduction 

Our goal is to study, from a convex geometric point of view, domino tilings of rectangular 
boards. They are in natural bijection with perfect matchings of grid graphs. Namely, we consider 
their associated perfect matching polytopes and present some of their properties. In particular, we 
characterize when these polytopes are Gorenstein. (We will define all these notions shortly). We 
also introduce the notion of domino stackings and present some results and several open questions. 

The m X n grid graph Q(m,n) is defined with vertex set G : < i < n, < j < m} 

such that two vertices and {i' are adjacent if and only if \i — i'\ + \ j — j'\ = 1. The m x n 
torus graph ^7^(m,n) consists of the same vertex and edge set as Q(m,n) with the additional 
edges {(0, j), (n — : < j < m} and {(i,0), {i,m — 1) : < i < n}. We use the convention 
that V{G) and E[G) denote the vertex set and edge set of a graph G. 

Given a graph G, M C E{G) is a perfect matching if every vertex of G is incident with exactly 
one edge of M. A generalization of a perfect matching is the notion of a magic labelling of sum 
i, which is a function E{G) Z>o such that for each vertex the sum of the labels of the edges 
incident to v equals t. Perfect matchings are magic labellings of sum t = 1. A natural question 
is how many magic labellings of sum t a given graph G has. In this paper, we are interested in 
the case when G = G{m,n), and denote the number of magic labellings of sum t by T{m,n,t). In 
particular, we can fix any two parameters, let the other vary to get a sequence, and encode this 
sequence in a generating function. 

The perfect matching polytope V associated to a graph G is defined to be the convex hull 
in R^^*^) of the incidence vectors of all perfect matchings of G. There is a natural identification 
between points in V and weighted graphs. We denote by V{ra, n) and VTim, n) the perfect matching 
polytopes of Q{m,n) and QT{m,n), respectively. 

For an integral polytope V (i.e., the vertices of V have only integer coordinates), let L'p{t) be the 
function that counts the number of integer points in tV := {tx : x G V} for t G Z>o. It was proved 
by Ehrhart [IT] that L'p{t) agrees with a polynomial with constant term 1 for all positive integers. 
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See also [4] for a modern treatment of Ehrhart theory. It is easy to see (as we will show) that the 
magic labellings of sum t of Q{m,n) correspond bijectively to the integer points of tV{m,n). It is 
a well-known result (see [2H Chapter 4.3]) that if p(t) is a polynomial, then Ylt>oP(^)'^^ evaluates 
to a proper rational function (i.e., the degree of the numerator is strictly less than the degree of 
the denominator) of z as follows. 

for a polynomial h{z) = h^z^ + hk-iz^~^ + • • • + /iq where k < d and hk ^ 0. We are interested 
in properties of the sequence (/iq, /ii, . . . , /ifc), which we call the Ehrhart /i-vector of V{m,n). 
For instance, Stanley showed (for general integral polytopes) that these numbers are nonnegative 
integers [20j . We characterize the values (m, n) for which V{m,n) is Gorenstein, that is, the 
Ehrhart /i-vector is palindromic {hj = hp^-j) [SIE]- equivalent formulation for the Gorenstein 
property is that there exists an integer k such that L-po(t) = L-p(t — k) for all t > A; and Lpo(i) = 
for t < k. Here L-poit) denotes the number of interior integer points of tV. In this case, we say 
that V is Gorenstein of index k. 

To warm up, we show in Section [2] that the Gorenstein property holds for VTifn,n) for some 
values. We can conclude this immediately from the following more general result. 

Proposition 1. Let G be a k-regular bipartite graph with ^{V{G)) even. Then the perfect matching 
polytope V of G is Gorenstein of index k. In particular, we have the functional identity Lpo(t) = 
L-p{t - k) for t > k. 

Corollary 2. Assume m <n and n > 2 is even. 

(1) 'Pt'(1)2) is a point. 

(2) PT(l,n) and Pj'(2, 2) are Gorenstein of index 2. 

(3) 7^7(2, ra) is Gorenstein of index 3. 

(4) If m > 2 is even, then VT{m^n) is Gorenstein of index 4- 

It is worth noting that Tagami [24j has shown that these and 'Pt{2, 3) and Vt{2, 5) are all of the 
perfect matching polytopes of torus graphs which are Gorenstein; the latter two are Gorenstein of 
index 4 and 3, respectively. 

For the non-torus grid graphs, we prove the following classification in Section [3l which is our 
first main result. 

Theorem 3. Assume m < n. The polytope V = V{m,n) is Gorenstein of index k if and only if 
one of the following holds: 

(1) m = 1 and n is even, in which case V is a point, 

(2) m = 2, in which case k = 2 if n = 2, and k = 3 for n > 2, 

(3) m = 3 and n is even, in which case k = 5, or 

(4) m = n = 4, in which case k = 4. 

Thus, in precisely these cases, V{m, n) has a palindromic Ehrhart /i-vector. We can say more 
(Section H]): 

Theorem 4. If Vim, n) is Gorenstein then V{m,n) has a unimodal Ehrhart h-vector. If m and n 
are both even, then VTim,n) also has a unimodal Ehrhart h-vector. 

(The sequence ho, /ii, . . . , /i^ is unimodal if there exists j such that /iq < • • • < ^j-i ^ ^ 
hj+i > • • • > hk-) Proposition [9] in Section [2] gives the dimension of V{m,n), so coupled with this 
result, we have an infinite list of Gorenstein polytopes. 
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The second theme of this paper concerns domino tihngs. A domino tiling of an mxn rectangular 
board is a configuration of dominos such that the entire board is covered and there are no overlaps. 
An excellent survey on general tilings can be found in [1]. There is an obvious correspondence 
between domino tilings and perfect matchings of Q{m,n). A natural question is whether there 
is an analogue of magic labellings of sum t of G{m,n). One possible answer is given by domino 
stackings. A domino stacking of height t of an m x n rectangular board is a collection of t 
domino tilings piled on top of one another. Every such domino stacking gives a magic labelling of 
sum t of Q{m,n) in a natural way. In Section [5] we establish a simple connection between magic 
labellings and domino stackings by showing that this natural map is surjective via some general 
results for lattice polytopes: 

Proposition 5. Every magic labelling of sum t of Q{m,,n) can he realized as a domino stacking of 
height t of an m x n rectangular hoard. 

We should remark that this result is a consequence of Hall's marriage theorem, but that it will 
also follow from properties of the perfect matching polytopes. 

It was shown by Klarner and Pollack [15j that for t = \ and fixed m (or, by symmetry, fixed 
n), the generating function X]n>o -^C"^' -'-)-^" ^ proper rational function in z. In Section [5l we 
slightly modify this proof to obtain the same result for general t. 

Theorem 6. For fixed m and t, the sequences (T{m,n,t))n>o o.'^d (T{m,n, l)*)n>o Q^'^e given by a 
linear homogeneous recurrence relation. 

Finally, in Section [U we give a geometric approach to domino tilings. Our main result in this 
section is a new, geometric proof of the following: 

Theorem 7 (Propp [16j). If T{m,n,l) counts the numher of ways to tile an m x n rectangular 
board with 2x1 dominos, then the sequence (T{m,n, l))n>o given by a linear recurrence relation, 
and furthermore, this recurrence relation satisfies the reciprocity relation 

T{m, n, 1) = (-1)" T(m, -n - 2, 1) 

if m = 2 mod 4, and 

T{m, n, 1) = T(m, -n - 2, 1) 

otherwise. 

A remark: we interpret negative arguments to a recurrence relation to mean that the formula for 
the recurrence relation is run backwards. It was shown by Propp that even if one uses nonminimal 
recurrence relations, these values at negative numbers are well defined. 

2. Basic Properties of perfect matching polytopes 

Let G be a graph. For x G R^^*^) and S C V{G), we denote by d{S) the set of edges that are 
incident to exactly one vertex in S, and x{d{S)) is the sum of the weights of the edges in d{S) in 
the weighted graph associated with x. The following theorem gives an inequality description for V. 

Theorem 8 (Edmonds [9j). Let G he a graph with an even numher of vertices. A point x = 
{xe : e G E{G)) lies in the perfect matching polytope of G if and only if 

(1) Xe > for all e G E{G), 

(2) x{d{v)) = 1 for all v G V{G), 

(3) x{dls)) > 1 for all S C V{G) of odd size. 

A complete characterization of graphs for which condition ([3]) is redundant is given in [8j. In 
particular, the condition is redundant for bipartite graphs. 
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Proof of PropositionUl As we just remarked, condition ^ in Theorem [8] is redundant for bipartite 
graphs. Thus L-pit) counts the integer points x = (xg : e G E{G)) that satisfy 

(1) Xe > for all e G E{G), 

(2) x\d{v)) = t for all v G 

whereas the interior points in V are those points satisfying (1) with strict inequality. In both cases, 
each equation in (2) involves exactly k variables, from which we deduce that L-po (t) = if 1 < t < /c 
and L'p°[t) = Lp{t — k) for t > k. □ 

An important piece of data associated with a polytope is its dimension. Fortunately, for grid 
graphs there is a simple formula for the dimension of its perfect matching polytope. The situation 
for torus graphs is complicated only in the sense that there are several cases to consider. 

Proposition 9. Suppose that mn is even. 

(1) dim P(m, n) = (m — 1) (n — 1) . 

(2) If n > 2 is even, then dim 7^7^(1,71) = 1. 

(3) If n > 2 is even, then dim'PT(2,n) = n + 1. 

(4) If n > 1 is odd, then dim7^T(2, n) = n. 

(5) If m > 2 and n > 2 are both even, then dim7^7^(m, n) = mn + 1. 

(6) If m > 2 is even and n > 1 is odd, then dim 'Pr(m,n) = mn. 

The proof of this statement follows easily from Theorem [10] below. But we first need some 
definitions. A graph G is said to be matching covered if every edge of G belongs to a perfect 
matching of G. A brick is a 3-connected and bicritical graph (i.e., removing any two vertices results 
in a connected graph that has a perfect matching). A method for decomposing general graphs into 
bricks is given in |10j . and the number of such bricks obtained is denoted by B{G). We shall not 
need this machinery; for our purposes it is enough to know that B{G) = if G is bipartite, and 
that B{G) = 1 if G is a brick. There is no overlap since the definition of a brick prevents it from 
being bipartite. 

Theorem 10 (Edmonds-Lovasz-Pulleyblank [TO]). Let G he a matching covered graph. If V is 
the perfect matching polytope of G, then dimP = #{E{G)) - #{V{G)) + 1 - B{G). 

Proof of Proposition Both Q{m, n) and Qxim, n) are matching covered if m > 1 and n> 1. The 
number of edges of Q{m, n) is m{n — 1) + n{m — 1), and the number of vertices is mn. In the case 
that either m = 1 or n = 1, V{m,n) consists of a single point, which has dimension and agrees 
with the formula. 

When m and n are both even, QT{m,n) is bipartite, so it is enough to count edges. For m = 2, 
the number of edges is 3n, and when m > 2, the number of edges is 2mn. 

In the other cases when m is even and n > 1 is odd, ^^(m, n) is not bipartite, but is 3-connected 
and bicritical if n > 1. So in this case, dim'Pj'(m, n) = ^[E{G)) — ^{V{G)), and we just need to 
count edges, which we've done above for n > 1. For the case m > 2 and n = 1, there are exactly 
two perfect matchings of Qj'{m, 1), which gives dim'P5"(m, 1) = 1. □ 

3. Characterization of Grid Graphs with Gorenstein Polytopes 

Now we give a lemma that will be used in the proof of Theorem [3l 

Lemma 11. Let n > 2 be even. Given a magic labelling of sum t of Q{3,n), let ai, hi, and Ci 
denote the values of the i^^ edge in the top, middle, and bottom rows, respectively. Then for i even, 
Ci = bi — a, and for i odd, Ci = t — Oi + bi. 
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Ci Ci+i ei_|_2 

Figure 1 . Proof of Lemma [TTJ 

Proof. We refer to Figure [TJ Here we are assuming that i is even and 1 < z < n — 4. Following the 
constraints given by the vertices, we get di+2 = t — ^i+i — cLi+2, which implies 

ei+2 = t — {t — aj+i — aj+2) — ^i+i — 
= Oj+i + — — ■ 
By induction, Cj+i = t — Oj+i + so 

Cj+2 = t — (t — aj+i + — (oj+i + aj+2 — — 64+2) 
= fei+2 — ai+2 • 

We play the same game for Ci+3. Solving some equations gives (ii+3 = t — aj+2 — flj+s and 

ei+3 = t — (t — aj+2 — Cj+s) — 6j+2 — 
= aj+2 + ai+3 — 6j+2 — 6i+3 • 
Using this identity, we finally conclude 

Cj+3 =t — Ci+2 — Cj+S 

= i — {bi+2 — aj+2) — (aj+2 + ai+3 — 6j+2 — fei+s) 
= t — aj+3 + 6i+3 . 

However, we have already taken care of the base case, as setting variables with negative index to 
be is equivalent to not having them there at all, and the proof proceeds in the same way. □ 

Now we turn to the proof of Theorem [3l Recall its statement: V = V{m, n) is Gorenstein if and 
only if one of the following holds: 

(1) m = 1 and n is even, 

(2) m = 2, 

(3) m = 3 and n is even, or 

(4) m = n = 4. 

Proof of Theorem\^ In order to show that a polytope V is Gorenstein of index A;, we prove two 
things. First, we show that kV is the smallest dilate of V to contain an interior integer point, and 
that this point is unique. Second, we show that if p is the unique integer point in kV° , then for 
any integer point x G tV° , Xi > Pi. If our graph is bipartite, then this is enough to conclude that 
X — p & {t — k)V using the hyperplane description given by Theorem [8l This gives an obvious 
bijection which implies that L-po(t) = — k). 

If m = 1 and n is even, then there is exactly one matching of Q{m, n). Then 7-" is a point, so the 
Ehrhart series is Ehr-p(z) = and the statement follows. 

We now consider the case m = 2. If n = 2, P is a 2-regular graph and Proposition [1] applies. 
For n > 2, note that the polytope tV has no interior points if t < 3. For t = 3, an interior point is 
given in Figure [2l The vertices of degree 3 force their adjacent edges to have weight 1, so it follows 
that this is the only such interior point. For t > 3, the edges of an interior point of tV must have 
weight > 2 if they have weight 2 in the figure. If not, then their adjacent edges have weight t — \ 
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Figure 2. The interior point of 3V for Q{2,n). 



^3^1^3^ ^3^1^3^ ^3^1^3^ 

211111111112 

211111111112 
1^3^1^34 1^3^1^3^ i^3^1^34 

Figure 3. The interior point of bV for G{3,n) and n even. 



and they are adjacent to vertices of degree 3, which gives a contradiction. So V is Gorenstein in 
this case. 

For the case m = 3, we may assume that n > 2. The existence of vertices of degree 4 means 
that tV has no interior points for t < 4. To get an interior point for t = 4, consider a vertex v of 
degree 2. If the edges adjacent to v have weight 1 and 3, then the edge of weight 3 is adjacent to a 
vertex of degree 3, which cannot happen. Thus, they must both have weight 2. However, there is 
a vertex w of degree 3 that is adjacent to two vertices of degree 2. Since the edges that they share 
both have weight 2, there is no vahd weight for the third edge adjacent to w, so also has no 
interior points. Using the notation of Lemma [TTl we get an interior point of bV with the following 
labelling: Oj = Cj = 3 if z is odd and = Cj = 1 otherwise, 5j = 1 if i is odd and 5j = 2 if i is 
even, (ij = ej = 2ifz = lorz = n + l and dj = = 1 otherwise. We illustrate this interior point 
in Figure [3l To see this is the only interior point of 5V, we consider an alternative labelling. No 
edge can have weight 4, so the weights for vertices of degree 2 must be 2 and 3. Note that the 3 
and 2 must be assigned as they are in Figure [31 Otherwise, the vertex of degree 3 on the far left 
will have sum greater than 5. We will show that for t > 5, Oj, Cj > 3 if i is odd and that 6j > 2 if 
i is even, which implies the uniqueness of the given interior point. By Lemma [TT| we know that 
Ci = bi — ai if i is even and Ci = t — ai + hi otherwise. For the first claim, suppose otherwise. Then 
3 > t — tti + bi > t — tti + l, which implies > t — 2, or Oj > t — 1. However, is adjacent to a vertex 
of degree 3 if i > 1, so this is a contradiction. For the second claim, if 6i = 1, then Cj = 1 — ai, 
which means Cj < 0, which cannot happen. Thus, we have shown uniqueness of the interior point 
in bV. But we have done more. From the inequalities shown, we get for free the functional identity 
L-po(t) = L'p{t — k), giving that V is Gorenstein. 

When m = n = 4, there are no interior points for tV when t < 3, and for t = 4, a unique interior 
point is given by labelling the eight edges incident to vertices of degree 2 with weight 2 and the 
rest of the edges weight 1. For a general interior point in tV when t > 4, the edges adjacent to the 
corner edges must have at least weight 2, so V is Gorenstein. 

For the other (m, n) of interest, we split them up into two cases. In the first case, m > 4 and 
n > 4 are both even, and in the second case, m > 4 is even and n > 4 is odd. 

First suppose that m > 4 and n > 4 are both even. There are no interior points in tV for 
t < 4. We give an interior point for AV as follows. Assign weight 2 to the edges {(0, i), + 1)}, 
{(m — 1, i), (m — 1, z + 1)}, {(i, 0), (i + 1, 0)}, and {(i, n — 1), (z + 1, n — 1)} when i is even, and assign 
1 to all other edges. We illustrate this for ^(6,6) in Figure HI This example is big enough to be 
instructive. To see this is unique, note that the weights of the edges of the corner vertices must be 
2, and the weights of the edges incident to a vertex of degree 4 must be 1. This, however, forces the 
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2 11112 
1^2^1^2^1^24 

Figure 4. The interior point of AT for C/(6, 6). 

|^2^2^1^2^2^ 

3 12 2 13 

12 112 1 



12 112 1 

3 12 2 13 

i—2 i 2 i 1 i 2 i 2— i 

Figure 5. An interior point of 5V for ^(6,6). 

rest of the values on the edges incident to vertices of degree 3. Now consider the fohowing interior 
point of bV. 

• Assign weight 2 to edges of the form {(i, 0), (i, 1)}, {i, and {{i, n— 2), {i,n—l)}, 
for i G {0, m — 1} and j odd. 

• Assign weight 3 to edges of the form 

• Assign weight 2 to edges of the form 

• Assign weight 2 to edges of the form 

• Finahy, assign weight 1 to ah other edges. 

We iUustrate this for ^(6,6) in Figure [5] (here (0,0) is the bottom left vertex and (5,5) is the top 
right vertex). In order for V to be Gorenstein, we must have L-po(t + 4) = L-p{t). Note that for 
any point of V, we may add the interior point of 47^ to get an interior point of bV. However, all 
points acquired through this must have weight > 2 for the edge {(0,2), (0,3)}, but we have given 
an interior point of 5V such that this edge has weight 1, which implies L-po(5) > L-p(l), so the 
functional identity is not satisfied. 

In the case that m > 4 is even and n is odd, tV does not have an interior point for t < 4. To see 
that this is true for t = 4, the edges incident to the corner vertices must have weight 2. This forces 
a labelling as in the case n is even, but it will not work because n is now odd. In other words, if we 
label the outside and work toward the middle, there will be no satisfactory values for the middle 
edges. We now give an interior point for bV. 

• Assign weight 3 to all edges of the form {{i, j), {i, j + 1)} where i G {0, n — 1} and j is even. 



(i + 1, j)} for i even and j G {0, n — 1}. 
(i + l,j)} for i odd and j G — 2}. 
(i + 1, j')} for i even and 2 < j < n — 3. 
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Figure 6. An interior point of bV for Q{6,7). 

• Assign weight 2 to edges of the form {i + 1, j)} where j € {0,m — 1} and i ^ 1. 

• Assign weight 2 to edges + 1)} where i € {1,2} and j is even. 

• Assign weight 2 to edges + 1)} where 3 < i < n — 2 and j is odd. 

• Finally, assign weight 1 to all other edges. 

This labelling is illustrated for Q{6, 7) in Figure [6] (again, (0, 0) is the bottom left vertex and (5, 6) 
is the top right). Since this graph labelling is not symmetric, we can get another interior point 
by "flipping" the graph (i.e., the weight for edges + 1)} and {(i,j), {i + 1, j)} are the 

weights for {{n — 1 — {n — 1 — i, j + 1)} and {(n — 2 — i,j), (n — 1 — i,j)}, respectively, from 
the interior point given previously). Thus, does not have a unique interior point so V is not 
Gorenstein. □ 

We give an example of the Ehrhart polynomial and Ehrhart series of the perfect matching 
polytope of a grid graph. We used polymake [l2] for conversion from the hyperplane description 
given in Theorem [8] to vertex descriptions, and normaliz [6] for the computation of the Ehrhart 
functions. 



Example. The Ehrhart functions of P(3, 4) are 

-^•P(3,4)(i) 



{t + l)(t + 2){t + 3){t + 4)(t2 + 5t + 5) 



Ehrp(3_4)(z) 



120 

+ 4z + 1 



{1-zy ■ 

The smallest example of a non- Gorenstein polytope arises from 7^(4,5), which has Ehrhart series 
21^*^ + 760z^ + 5919z^ + 15578^^ + 16432^"^ + 73562^ + 1339^^ + 82z + 1 



Ehr 



•P(4,5) 



(1 



U3 



4. Unimodular Triangulations and Unimodality 



We remind the reader that a triangulation of a polytope P is a set of simplices A whose union 
is V such that the intersection of any two simplices is a face of both. An integral simplex with 
vertices vq, . . . ,Vn is unimodular if {vi — vq, . . . ,Vn — vq} forms a basis of Z" as a free Abelian 
group, and a triangulation is unimodular if each simplex is unimodular. 

One might ask if V{m, n) has a unimodular triangulation. Fortunately, the answer is yes. Given 
an ordering r = (vi, . . . ,Vn) of the vertices of a polytope V, the reverse lexicographic trian- 
gulation is defined as follows. For a single vertex v, A^ = {v}. In general, consider all facets 
(codimension 1 faces) of V that do not contain w„. For each such facet F, there is an ordering 
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of the vertices of F induced by r that gives a reverse lexicographic triangulation of F. For each 
simplex in this triangulation, take the convex hull of its union with w„. The union of all such convex 
hulls is A^. A polytope is said to be compressed if the reverse lexicographic triangulation with 
respect to any ordering of its vertices is unimodular. The following result immediately guarantees 
the existence of unimodular triangulations. 

Theorem 12 (Santos [18j, Ohsugi-Hibi jilSj, Sullivant ^23j). If V is a 0/1 polytope (i.e., each 
coordinate of each vertex is either or 1) defined by the linear equalities and inequalities 

A X = u and < Xj < 1 , 

where A is an integer-valued matrix and u is an integer-valued column vector, then V is compressed. 

As a bonus, we obtain Theorem [H that is, for the {m,n) discussed in Theorem [3l the Ehrhart 
/i-vector is unimodal. This is an immediate consequence of Theorem [12] combined with the following 
result. 

Theorem 13 (Athanasiadis-Bruns-Romer [21 [7]). If V is compressed and Gorenstein, then the 
Ehrhart h-vector ofV is unimodal. 

5. Domino Stackings and Magic Labellings 

There is a natural many-to-one function from domino stackings of height f of an m x n rectangular 
board to the magic labellings of sum t of Q{m, n). It is conceivable that some magic labelling cannot 
be realized as a physical stacking. However, we have shown earlier that V{m, n) has a unimodular 
triangulation, which implies that it is a normal polytope. (An integral polytope V is normal if 
every integer point in tV can be written as a sum of t integer points in V.) In particular, this 
means that every magic labelling of sum t of Q{m,n) must be the sum of t perfect matchings, so 
our function is a surjection, and this is what we claimed in Proposition [5j 

Recall that for fixed m, the generating function X^„>o T(m, n, evaluates to a proper rational 
function in z, which means that the sequence {T{m, n, l))ri>o is given by a linear recurrence relation. 
Running these recurrence relations backwards, one obtains values for T(m, n, 1) when n < 0. It 
was shown in [T6] that these values for n < are the same even if the recurrence relation used is 
not minimal. A reciprocity relation was also given, which is the content of Theorem [71 

This is equivalent to a result in [22] which gives properties of the numerator and denominator 
polynomials in the rational function ^^>g T(m, n, 1)^*. We would like an analogous result for 
magic labellings of sum t and domino stackings of height t. Recall that Theorem [6l states that for 
general t, if one fixes m, then there is a recurrence relation in n for the function T(m,n,t). The 
proof that we give is essentially the same as the one found in [15j with only slight modifications. 

Proof of Theorem\^ Given a magic labelling of sum t for Q{m,n), consider its canonical planar 
drawing. We can uniquely encode it with n + 1 vectors cq, . . . , of length m in the following way. 
For 1 < i < n — 1, the j^^ entry of Cj is the label of the j^^ horizontal edge counting from the top of 
the graph in the i*^ column starting from the left and going right. We define cq and c„ to be zero 
vectors of length m. 

We now form a directed graph whose vertex set is {0, 1, . . . ,4}*". For two vectors u,v & 
{0, we form the edge {u,v) if and only if there is some n such that there is a magic 

labelling of sum t for Q{m, n) encoded by the vectors cq, . . . , c„ with Ci = u and q+i = v for some i. 
Given this, there is a natural bijection between the magic labellings of sum t of G{rn, n) and walks 
of length n in from {0}"* to itself. 

Order the vertices of by interpreting each vector as a number in base t -\- 1, and let be 
the adjacency matrix of S^m- The number of walks of length n from {0}™" to itself is the (0, 0) entry 
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of A^. By the Cayley-Hamilton theorem, there is a polynomial relation 

Ai + ftd-i^^-^ + • • • + biArn + bo = 

for d = {t + 1)™" and integers 6j. Multiplying both sides by A!^'^ and looking at the (0,0) entry in 
each term gives a recurrence relation for T{m,n,t) involving the coefficients 6j. 

In the case of domino stackings of height t, we can get a similar result. Instead of vectors in 
{0, . . . , t}^ making up the vertices of "^rm we encode each domino stacking as a kxt zero-one matrix 
and play the same game to get a recurrence relation for the sequence (T{m,n, l)*),i>o. □ 



As an immediate corollary, we now know that for fixed m and t, the generating function 

P(z 



^„>Q T(m, n, evaluates to a rational function ^fy- The next step is to understand what 



properties P and Q have. 

The sequence (T(2,n, l))n>o is a shift of the Fibonacci numbers, so this gives a nice proof that 
powers (for a fixed exponent) of Fibonacci numbers can be encoded by a linear homogeneous 
recurrence relation. This was previously done in [T7], which also gives recurrence relations for the 
coefficients of the recurrence relations. 

6. Convex Geometric Considerations of Domino Tiling Reciprocity 

In this section, we give a convex geometric proof of why reciprocity exists for the number of 
domino tilings of the m x n rectangular board for fixed m using the following theorem which gives 
a closed form. We make some remarks on reciprocity of domino stackings of fixed height t and 
width ni at the end of the proof. 

Theorem 14 (Kasteleyn [H]). The number of domino tilings of the m x n rectangular board is 
given by 

lm/2] ^n+l _ -n+l 

T(m,n,l)= 4 — 

j=i ~ 

where 



C, = COS h A / 1 + cos^ 



Cj = COS ■ 



m + 1 V m + 1 

7 ~ A/ 1 + cos"= . 

m + 1 V m + 1 

We remark here that the above theorem is in a different form from the one found in [14j. The 
original formulation gives a product from j = 1 to j = [m/2j , and designates r(m, n, 1) = when 
mn is odd. However, when mn is odd, c""*"^ — c""*"^ = for j = \m/2\ , and for m odd and n even, 
it is Cj — Cj, so our formulation is equivalent. 

First note that CjCj = — 1 for all j. The key observation is that 

^n+l _ -n+l 

=c» + cr^c, + ... + c,c^-^ + c-, 

and that the sum on the right is an evaluation of a certain multivariate generating function ([*]) 
below. That is, consider the polytope V = conv{(l, 0), (0, 1)}. Let i := \m/2]; we are interested 
in the ^-fold Cartesian product of V with itself. Let C C R^^+^ be the cone over V. To be precise, 
C = {{x, A) G X R : A > 0, x € XP}. We consider the integer point transform of C, which is 
the following formal Laurent series in 2^ + 1 variables. 

(*) c^c(x) := Yl ^ 

yGCnZ2'^+i 
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where := ■x^^'x}^ ■ ■ ■ x^2l+i- -'-'^ case of a simplicial cone C, if yi, . . . , are the generators of 
C, define the fundamental parallelepiped 

n(C') := {Aiyi + • • • + XkVk : < Ai < 1} . 

Then the integer point transform (Tc'(x) is a rational function of the form (see, e.g., [U Chapter 3]) 

c^n(C')(x) 



<7c'(x) 



(1 - xJ^i) • • • (1 -xffe) 



We can triangulate V (using no new vertices), and this triangulation induces a triangulation of C 
into simplicial cones. The generators of each will be of the form XjjXjj ■ ■ ■ a^j^ where {ii, • • • > h} ^ 
{1, . . . , 2f\ contains exactly one of 2i — 1 and 2i for each i = Conversely, each such subset 

will appear as the generator for at least one of the simplicial cones in the triangulation of C. 
By adding the rational functions for these simplicial cones and playing inclusion-exclusion on the 
intersections, the integer point transform of C is 

£^c(x) - 



where the product runs over all ^-subsets of {1, ... , 2C\ that contain exactly one of 2i — 1 and 2i for 
each i = and P is some multivariate polynomial. We make the substitutions X2i~i i— > Q, 

X2i ^ Ci, and X2i+i i— > z. The polytope nV can be described as the set of nonnegative solutions to 
the equations xi + • • • + X2i = in and X2i-i + X2i = n for each z = 1, . . . , Thus, using Theorem[li 

I 

Cr„p(ci, Cl, . . . , Q, Q, Z) = JJ(C^ + C^-^Cj + • ■ ■ + CjC^-^ + C^) 

i=i 

= r(m, 71,1)2". 

From this substitution, 

(Tc(ci,ci, . . . ,Q,Q,z) = ^r(m,n, l)z"' 

n>0 

P(ci,ci, . . . ,Q,Q,2;) 



n,(i-di^4^---d^) 



where / runs over all possible binary strings of length £, and ij denotes the j^^ bit of I, and (f- = Cj 
and = Cj. 

To proceed, we use the following theorem due to Stanley which gives a relation between the 
integer point transform of a rational cone and the integer point transform of its interior. 

Theorem 15 (Stanley [19]). IfC is a rational cone in R"^ whose apex is the origin, then as rational 
functions, 



ac{z^\...,z-,^) = {-lt^^aco{z^,...,Zd), 



where C° denotes the relative interior of C. 
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Equipped with this fact, we have 

CJc°(ci, Ci, . . . , Q, Q, z) = (-1)^+Vc(-Ci, -Ci, . . . , -C£, -Q, Z"^) 

= (-l)^+Vc(ci,Cl,...,Q,Q,(-lfz-l) 

(1) = (-1)^+1 ^(-l/"r(m,-n,l)z« 

n<0 

= (-l)'^^(-l)^"r(m,-n,l)z". 

n>0 

The last equahty follows from the identity (of rational functions) 

^/(n)z" + ^/(n)z" = 0, 

n>0 n<0 

where f{n) is given by a linear recurrence relation. Finally, we use the fact that C is a Gorenstein 
cone (i.e., adding some vector to the cone allows us to enumerate the integer points in its interior). 
This means that 

o-c°(ci,ci, . . . ,ce,Ci,z) = c\c\ ■ ■ ■ cecez'^ac{ci,ci, . . . ,ce,ce,z) 

(2) = (-l)V^T(m,n,l)z". 

n>0 

Putting both sets of equalities together and matching coefficients gives 

r(m, n, 1) = (-1)^" r(m, -n - 2, 1) . 
To recover the result in [16], note that when m = 3 (mod 4) or m = (mod 4), ^ is even, so 

T(m, n, 1) = T(m, —n — 2,1). 

In the other cases, i is odd, so 

T{m, n, 1) = (-1)" T{m, -n - 2, 1) . 

In particular, for m = 1 (mod 4), T(m, n, 1) = when n is odd, so this is equivalent to 

T(m, n, 1) = T(m, —n — 2,1), m = 1 (mod 4) , 

which gives the desired reciprocity relations. 

For general t, raise both sides of the reciprocity equations to the t^^ power to get reciprocity 
relations for domino stackings of height t. Alternatively, we could run through the above proof 
with the following small modifications. Consider the t-fold product of V with itself. The cone over 
is now {ti + l)-dimensional, so the last term in equation ([1]) becomes 

(-l)«^(-l)«"r(m,-n,l)V, 

n>0 

and the last term in equation ([2]) becomes 

(-l)"z2^r(m,n,l)*z". 

n>0 

Combining these two gives the general reciprocity relation 

T(m, n, 1)* = (-1)"" r(m, -n - 2, 1)*. 
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7. Further Directions 

A natural question to ask is if unimodality of the Ehrhart /i- vector extends to the (m, n) such 
that V{m,n) is not Gorenstein. Unfortunately, computing examples larger than ^(4,5) becomes 
difficult, so there is not much evidence supporting whether they are or are not unimodal. 

We defined Q{m,n) having vertex set €Z^:0<i<m, 0<j<n} with {(«, j), 

an edge if and only if |i — i'| + \j — j'\ = 1, but there is no reason to restrict to 2-dimensional lattices. 
We can define /c-dimensional grid graphs Q{mi,m2, • • • , rti]^) having vertex set {(xi, . . . , rr^) E Z'^ : 
< Xj < nii} such that {(xj), {vi)} is an edge if and only if Yl\=i \xi~yi\ = 1- What further things 
can one say about these higher dimensional analogues? It would make sense that very few of these 
graphs have Gorenstein polytopes. The perfect matching polytope of the cube graph Q{2, 2, . . . , 2) 
is Gorenstein for arbitrary k, since this is a /c-regular graph. Some questions one can ask is if 
there is a largest k such that ^(2, 2, . . . , 2) is the only graph with a Gorenstein polytope (assuming 
that TTij > 1). If not, is there a largest k such that the number of nontrivial Gorenstein polytopes 
V{mi, . . . ,mk) is infinite? We can also ask about the function T{mi, . . . ,mk;t) that counts the 
number of magic labellings of sum t of Gimi, . . . ,mk). Since multi-dimensional grid graphs are 
bipartite, X]t>o-^("^i' • • • i''^k',t)z^ evaluates to a rational function. It is not difficult to see using 
the proof of Theorem [6] that J2m^>o ^^i^ii ■ ■ ■ > mk;t)z'^^ is also a rational function. 

One could also define two domino stackings to be equivalent if the layers of tilings of one is a 
permutation of the layers of tilings of the other. Is there a recurrence relation for this new set of 
stackings? New techniques would be needed as the adjacency matrix argument no longer works. 

As mentioned earlier, Riordan gives recurrence relations for coefficients of recurrence relations 
for the sequences given by powers of Fibonacci numbers in [17]. It is conceivable that the same 
phenomenon occurs for powers of T(m, n, 1) for general m. Different techniques from what is 
presented here would most likely have to be employed. The main problem arises in that the proof 
of Theorem [6] does not give minimal recurrence relations. For instance, in the case m = 2 and t = 1, 
which is a shift of the Fibonacci sequence, a 5-term recurrence relation is given by the proof, but 
only a 3-term recurrence relation is needed. This corresponds to the fact that if this non-minimal 
recurrence relation appears in the denominator of the rational function ^^-^q T(m, n, 1)*2;"', then 
it will not be in reduced terms. A further problem is that the adjacency matrix argument given 
above presents a recurrence relation that works for all of the entries of the matrix, but we are only 
interested in one specific entry. 

Along those same lines, for fixed t, the sequence T(m, n, t) for counting magic labellings of 
G{m,n) is given by a linear homogeneous recurrence relation if either m or n is fixed. It would be 
interesting to find a recurrence relation when m and n are both allowed to vary, and also to see 
if the coefficients that arise from these recurrence relations can also be shown to satisfy their own 
recurrence relation. 
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